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ABSTRACT 

The evolution of magnetic fields is studied using simulations of forced helical turbu- 
lence with strong imposed shear. After some initial exponential growth, the magnetic 
field develops a large scale travelling wave pattern. The resulting field structure pos- 
sesses magnetic helicity, which is conserved in a periodic box by the ideal MHD equa- 
tions and can hence only change on a resistive time scale. This constrains strongly 
the growth time of the large scale magnetic field, but less strongly the length of the 
cycle period. Comparing with the case without shear, the time scale for large scale 
field amplification is shortened by a factor Q, which depends on the relative impor- 
tance of shear and helical turbulence, and which controls also the ratio of toroidal 
to poloidal field. The results of the simulations can be reproduced qualitatively and 
quantitatively with a mean-field afl dynamo model with alpha-effect and the tur- 
bulent magnetic diffusivity coefRcients that are less strongly quenched than in the 
corresponding a^-dynamo. 



1 INTRODUCTION 



In astrophysical bodies such as stars and galaxies there is 
a strong magnetic field of large scale. Such fields have usu- 
ally significant magnetic helicity (e.g., Pevtsov et al. 1995, 
Berger & Ruzmaikin 2000). This is non-trivial, because mag- 
netic helicity is a conserved quantity and can only change if 
there is a flux of helicity through the boundaries, or through 
resistive effects which are however very slow. Although this 
has been known for some time, it is only recently that this 
has been identified as the fundamental reason for 'catas- 
trophic' quenching of the a-effect in mean-field dynamo the- 
ory (Blackman & Field 2000, Kleeorin et al. 2000). Simula- 
tions of non-mirror symmetric turbulence, which is proto- 
typical of flows producing a^-dynamos, have shown that a 
large scale helical magnetic fleld can only grow to its final 
(super-) equipartition field strength on a resistive time scale 
(Brandenburg 2001, hereafter referred to as B2001). 

One may be tempted to sweep the problem of helicity 
conservation under the carpet, because it has mainly been 
discussed in connection with rather idealised models. We be- 
lieve however that the problem is serious and quite general. 
In fact, it also applies to convection-driven dynamos and 
even to the case where the dynamo-generating flow is the re- 
sult of magnetic instabilities, as was found to be the case in 
simulations of accretion discs with a dynamo-generated large 
scale field (Brandenburg et al. 1995). This may be particu- 
larly surprising in view of the rather plausible expectation 
that the a-effect and turbulent diffusivity should be 'anti- 
quenched' and increase with increasing field strength (Easier 
et al. 1995, Brandenburg et al. 1998). If such a mechanism is 
to be successful, it must still obey helicity conservation and 



can hence only produce a field with vanishing net magnetic 
helicity. 

There is strong observational evidence that the solar 
magnetic field is indeed helical (Seehafer 1990, Pevtsov et 
al. 1995). These observations suggest negative current he- 
licity of the small scale fields in the northern hemisphere. 
Using a relation by Keinigs (1983), this implies a positive 
a-effect (Seehafer 1996), which is consistent with B2001. In 
order to produce finite net-helicity one must get rid of fields 
with opposite sign of magnetic helicity, either through dis- 
sipation (which is slow) or through selective losses through 
open boundaries. So far there is no evidence from simula- 
tions however that such losses involve fields of significant 
strength and opposite sign of magnetic helicity relative to 
those that remain in the dynamo-active domain (Branden- 
burg & Dobler 2001). 

The dynamo simulations that allowed addressing the 
question of the helicity constraint were all of a^-type, so 
there was no additional field amplification by shear. Thus, 
an outstanding question is therefore whether or not the he- 
licity constraint also plays a role in the presence of shear 
through which strong toroidal magnetic fields can be gener- 
ated without affecting the magnetic helicity. 

There are a number of working dynamos which have 
both open boundaries and shear (e.g., Glatzmaier & Roberts 
1995, Brandenburg et al. 1995), but those models are 
rather complex and use subgrid scale modelling, so one can- 
not straightforwardly define an effective magnetic Reynolds 
number. This makes a reliable assessment of the effects of he- 
licity conservation difficult. Nevertheless, it clearly remains 
one of the next important tasks to reconsider these or similar 
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simulations in the light of helicity conservation. In order to 
determine the relative importance of the various possibilities 
for relaxing the helicity constraint (shear, open boundaries, 
etc.) it is useful to consider each possibility in isolation. As 
a straightforward extension of the model of B2001 we con- 
sider here the inclusion of large scale sinusoidal shear, which 
allows us to retain the assumption of periodic boundary con- 
ditions. 

We have mentioned already that shear could be impor- 
tant for relaxing the helicity constraint because the toroidal 
field generated by stretching does not need to be helical and 
would hence not be subject to the helicity constraint. On the 
other hand, shear alone is insufficient for dynamo action: one 
needs an additional effect that regenerates poloidal (cross- 
stream) field from toroidal field (e.g. Moffatt 1978, Krause 
& Radler 1980). The main point of the present paper is to 
show that, even though much of the magnetic field ampli- 
fication is due to shear, which causes the field to be only 
weakly helical, the magnetic field is still subject to a (mod- 
ified) helicity constraint. More specifically, we shall show 
that it is no longer the large scale field as such which grows 
resistively, but rather the geometrical mean of the magni- 
tudes of the poloidal and toroidal mean fields. The reason is 
simple: large scale helicity measures essentially the linkage 
of poloidal and toroidal fields and must therefore be propor- 
tional to the product of the two. The constraint that helicity 
can change only on a resistive time scale can then be allevi- 
ated somewhat. This is because, now, for the same magnetic 
helicity, stronger toroidal fields are possible at the expense of 
weaker poloidal fields. Or conversely, equipartition strength 
large scale fields can be attained in times shorter by the ratio 
of toroidal to poloidal field strength. 



2 THE MODEL 

As in B2001 we adopt the MHD equations for an isothermal 
compressible gas, driven by a given body force /, which 
represents both shear and small scale driving; 



fk- 



k X {k X e) — i|fc|(fe x e) 
2fcVl - (fe ■ e)Vfc2 



(7) 



D Inp 



= -V ■ u, 



(1) 



Dw JxB /i ,„2 X ,. 

— = -CsVlnp-h + t^{Vu+^VV-u) + f, 2 

Dt p p 



OA 

dt 



u X B — rjfioJ, 



(3) 



where D/Dt — d/dt + u ■ ^ is the advective derivative, u 
is the velocity, p is the density, B — ^ x A is the magnetic 
field, A is its vector potential, J = V x B / no is the current 
density, rj is the magnetic diffusivity, and p, the dynamical 
viscosity. We adopt a forcing function / of the form 

f f turh f shear ^ (4) 

where 

fshcai = C'shcar — y siu X (5) 
P 

balances the viscous stress once a sinusoidal shear flow has 
been established, and 



/turb = Re{Af/j^^jj exp[ife(t) ■ x + i(t>{t)]}, 
is the small scale helical forcing with 



(6) 



Here e is an arbitrary unit vector needed in order to generate 
a vector k x e that is perpendicular to k, (j>{t) is a random 
phase, and = foCs(kCs/Sty^'^ , where fo is a nondimen- 
sional factor, k — \k\, and St is the length of the time step. 
As in B2001 we focus on the case where |fc| is around kf = 5, 
and select at each time step randomly one of the 350 possible 
vectors in 4.5 < |fe| < 5.5. 

We use nondimensional units where Cs = fci = po = 
^0 = 1- Here, Cs is the sound speed, fci is the smallest 
wavenumber in the box (so its size is 2n), po is the mean 
density (which is conserved), and po is the vacuum perme- 
ability. 

We are interested in the case where shear is strong com- 
pared with the turbulence, but still subsonic. In B2001 we 
used fo — 0.1 and found that the resulting Mach number of 
the turbulence was between 0.1 and 0.3, which is already so 
close to unity that there would be no room to accommodate 
sufficiently large shear which is still subsonic. Thus, we now 
choose fo to be ten times smaller, and we take fo — 0.01. 
During the saturated phase of the dynamo the resulting rms 
velocities in the meridional (xz) plane are now around 0.015. 
For the shear parameter we choose Cshcar ~ 1, which leads 
to toroidal rms velocities of around 0.6, which is about 40 
times stronger than the velocities in the meridional plane. 
The rms velocity from wavenumbers fc > 2 is 0.035, and 
this is also the value that we shall use for our estimates of 
the magnetic Reynolds number and the equipartition field 
strength. 

We choose a magnetic Prandtl number of 10, i.e. 
f-fiPov) ~ lOi f'lid use 77 = 5 X 10"*, so the magnetic 
Reynolds numbers based on the box size (= 2tt) is about 
400. The magnetic Reynolds number based on the forcing 
scale is about 80. The kinetic Reynolds number based on 
the forcing scale is only 8, so one cannot expect a proper in- 
ertial range. The turnover time based on the forcing scale is 
r = 40. In the following we denote by poloidal and toroidal 
components those in the xz-plane and the y-direction, re- 
spectively. 

As usual for these type of simulations with helical forc- 
ing, there is strong dynamo action at small scales ampli- 
fying an initially weak random seed magnetic field expo- 
nentially (on a dynamical time scale) to equipartition with 
kinetic energy. The poloidal field, which is strongly domi- 
nated by small scales, saturates early on (at t ~ 1000) at a 
level of about 0.010-0.015. The toroidal field saturates later 
(at t « 2000) at a level of about 0.2-0.3, and is then already 
dominated by large scales. 

We begin by discussing the resulting field structure at 
late times, turn then to the question of resistively limited 
growth of the large scale field, and finally make comparisons 
with aQ dynamo theory. 



3 FIELD STRUCTURE 

In Fig. we show images of the three field components in the 
meridional plane. Note that the toroidal field shows much 
smoother and larger scale structures than the meridional 
field components. Moreover, the toroidal field shows almost 
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Figure 1. Images of the three components of B in an arbitrarily 
chosen xz plane (first row), compared with the i;-averaged fields 
(second row) and the fourier-filtered y-averaged fields with |fc| < 
2, indicated by the subscript 'filt' (third row). 120"^ meshpoints, 
t = 6000. 



no variation along the y-directions: the toroidal average, By, 
(second row), is very similar to an individual meridional 
cross-section of By, first row. However, in contrast to the 
case without shear, where the mean fields showed system- 
atic variations only in one of the three coordinate directions 
(B2001), here the toroidal field varies with both x and z, 
consisting of a superposition of modes with kx = 1 and 
= 1. 

The toroidal component of the mean field displays dy- 
namo waves travelling in opposite directions at different 
x-positions, depending on the local sign of the shear. For 
X = —TV the local shear is negative and the dynamo wave 
travels in the positive z-direction, whilst for a: = the lo- 
cal shear is positive and the wave travels in the negative 
^-direction (at least after t = 4000); see Fig. ^. This is con- 
sistent with what is predicted from mean-field aQ dynamo 
theory (e.g. Yoshimura 1975) . The dynamo wave at a; = — tt 
is quite well established at t = 2000, but the behaviour at 
a; = is more complicated and a clear dynamo wave de- 
velops only after t — 4000. The cycle period at a; = is 
also longer than at x = —tt. This somewhat complicated 
behaviour suggests that the turbulence properties may not 
be homogeneous in x, which could be a consequence of the 
magnetic feedback. 

There is a systematic phase shift and a well-defined am- 
plitude ratio between By and B^; see Fig. ^. Note also that 
the dynamo wave is markedly non-harmonic. These are clear 
properties that can be compared with mean-field model cal- 
culations (Sect. 

Before we turn to the saturation of the field at the scale 
of the box we first want to assess the relative importance of 
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Figure 2. Space-time diagram of the mean toroidal field at 
X = —IT (negative local shear) and x = (positive local shear). 
Dark (light) shadings refer to negative (positive) values. Note the 
presence of dynamo waves travelling in the positive (negative) 
2-direction for negative (positive) local shear. 
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Figure 3. Evolution of Bx and By at x = 
that Bx has been scaled by a factor —100. 
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the different Fourier modes at different times. Thus, we plot 
in Fig. ^ the evolution of the power, \Bi{kj)\'^, in a few se- 
lected modes. Note that after t = 1700, most of the power is 
in the mode \By{k^)\^ , i.e. the toroidal field component with 
variation in the z-direction. Between t = 1700 until t ~ 3500 
the ratio of toroidal to poloidal field energies is around 10*, 
so Btor/ Bpoi ~ 50. At later times this ratio diminishes some- 
what. This may suggest that there is a growing contribution 
from a^-type dynamo action. This is also supported by the 
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Figure 4. Evolution of the power, \Bi{kj)\^, of a few selected 
Fourier modes. After t = 1700, most of the power is in the mode 
i.e. in the toroidal field component with variation in 
the z-direction. 




apparently independent evolution of the oscillatory fcz-mode 
and the non-oscillatory fc^-mode; see Fig. ^. 

In Fig. P we show two-dimensional power spectra of 
the three components of the mean field, B. (Here and else- 
where we denote y-averaged fields by a bar whilst angu- 
lar brackets are used for full volume averages.) Note that a 
strong toroidal field builds up first, and at later times the 
poloidal field components also gain significant power at the 
largest scale (i.e. at < 2). One should bear in mind, how- 
ever, that these spectra are for the mean fields. The three- 
dimensional power spectra of the non-averaged fields reveal 
that the poloidal fields are 'noisy' and possess significant 
power at the forcing wavenumber, k{- see Fig. ^. 

The small scale contributions to the poloidal field result 
from variations in the toroidal direction, as can be seen in a 
longitudinal cross-section; see Fig. ^, where we show images 
of the three field components in the yz plane. The figure 
shows that whilst the toroidal field is relatively coherent 
in the toroidal direction, the poloidal field components are 
much less coherent and show significant fluctuations in the 
y-direction. 

We now turn to the temporal evolution of the resulting 
large scale magnetic field that gradually emerges during this 
simulation. We begin by briefly reviewing the main results 
in the absence of shear (B2001). 
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Figure 5. Two-dimensional power spectra of the three com- 
ponents of the mean field. By (solid for the y component, and 
brolsen lines for the x and z components. The slope is given 
for comparison. 



4 RESISTIVELY LIMITED GROWTH ON 
LARGE SCALES 

In an unbounded or periodic system the magnetic helicity, 
{A ■ B), can only change if there is microscopic magnetic 
diffusion, rj, and finite current helicity, {J ■ B), 



{A-B) = -2r]{J ■ B). 



(8) 



dt 

In B2001 a possible configuration for the large scale mag- 
netic field was 

(COs(fcl2: + ip^) \ 
sin(fci2: + i^y) j , 
/ 



(9) 



which corresponds to a force-free magnetic field that varies 
in the z-direction, although variations in one of the other 
two coordinate directions, and with arbitrary phase shifts 
ifix (~ ipy), were also possible (B2001). -Bo ~ {B )^^^ is the 
amplitude, whose time dependence was found to be subject 
to the helicity constraint (B2001). 

The present case is different because of shear which 
tends to increase the toroidal field, but not the poloidal field. 
We model this by writing 



CBpol COs(fcl2 + ifx) 
Btor sm{kiz + ify) 




(10) 
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Figure 6. Three-dimensional power spectrum of the three field 
components. 120^ meshpoints, t = 6000. 
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{J-B}/ki ^TBtorBpoi~ki{A-B). 



(11) 



( J • B) = (J ■ B) + {j ■b)^0. 
This yields, analogously to B2001 
d 



(12) 



dt 



(Btor Bpol) = +2rjkt (Btor Bpoi) - 2riki\{j ■ (13) 



which yields the solution 

R R _ ^ r2 [i -2rik^(t-ts) 



(14) 



Figure 7. Images of the three components of B in an arbitrarily 
chosen yz plane. Note that Bx and show strong variations in 
t = 6000. 



where Bpoi and Btor are the amplitudes of the poloidal and 
toroidal field components. In addition to the z-dependence 
there can also be an a;-dependence of the mean field, which 
is natural due to the a;-dependence of the imposed shear pro- 
file. However, for the following argument all we need is the 
fact that the magnetic and current helicities are proportional 
to the product of poloidal and toroidal field magnitudes, 



The upper sign applies to the present case where the kinetic 
helicity is positive (representative of the southern hemi- 
sphere), and the approximation becomes exact if the field 
is indeed represented by Eq. (p^. 

Following B2001, in the steady case {A ■ B) — const, 
see Eq. (^), and so the r.h.s. of Eq. (^) must vanish, i.e. 
(J ■ B) = 0, which can only be consistent with Eq. (|ll]) 
if there is a small scale component, {j ■ b), whose sign is 
opposite to that of {J ■ B). Hence we write 



Figure 8. Growth of poloidal and toroidal magnetic fields on 
a logarithmic scale (upper panel), and product of poloidal and 
toroidal magnetic fields on a linear scale. For the fit we have used 



kf = 2, Ba. 



: 0.035, and eo = 1-3. 



where eo ~ \ {j 'l>)\/ {kiB'^^) is a prefactor, B^q is the equipar- 
tition field strength with B^q — /lo(pw^), and tn is the 
time when the small scale field has saturated which is when 
Eq. ( |l3| ) becomes applicable. All this is equivalent to B2001, 

except that {B ) is now replaced by the product Btor Bpoi. 
The significance of this expression is that large toroidal fields 
are now possible if the poloidal field is weak. 

In order to compare with the simulation we now define 



BtoT 



pol 



{Bx + B,y 



(15) 



Note that this definition generalizes that given in Eq. (|lO|). 
In Fig. 1^ we show the evolution of Stor and Bpoi and compare 
the evolution of the product Btor Bpoi with Eq. (|lj) . There 
are diflterent stages; for 1200 < t < 2200 and 3000 < t < 
3700 the effective value of kf is kf = 2 (because there are 
contributions from k^ = 1 and fcz = 1; see Fig. whilst at 
other times (2500 < t < 2800 and t > 4000) the contribution 
from kx = l (for 2500 < t < 2800) or fc^ = 1 (for t > 4000) 
has become subdominant and we have effectively kf = 1. 
This is consistent with the change of field structure discussed 
in the previous section: for 2000 < t < 3000 and around 
t — 4000 the By{kx = 1) mode is less powerful than the 
By{kz = 1) mode. 

We may conclude that the effect of the helicity con- 
straint is clearly identified in the present simulations. This 
is substantiated by the fit shown in Fig. ^ The episodes dur- 
ing which the field amplitude is below that obtained from the 
helicity constraint can be explained by temporary changes 
in the field geometry. 



5 ASTROPHYSICAL IMPLICATIONS 

The main result of this paper is a quantitative modification 
of the helicity constraint for dynamos in the presence of 
shear. With shear included the estimate for {B ) of B2001 
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is now to be replaced by the product Btor Bpoi ~ {B ) /Q, 
where Q = Btor/Bpoi 2> 1 and so (i? 
times, the exponential function in Eq 



(0 



U)Q. 



=t Bt,. For 
can be expanded: 

(16) 



In the case of efficient large scale dynamo action the small 
scale current helicity is very nearly equal to the normalized 
kinetic helicity, pqI^uj-u) (see also Brandenburg & Subrama- 
nian 2000), which in turn is approximately k{{pu'^). Since 
eo = (j • b)/{kiB'^^), this leads to eo ~ fcf/fci, which is 5 
in the present case. The value of eo that fits the simulation 
results best is only 1.3 (see Fig. ^), so the dynamo seems to 
be not fully efficient. This reduced efficiency could partly be 
explained by the fact that the actual field is not sinusoidal. 



as assumed in Eq. ( 
poloidal and toroida 



and that the phase shift between 
fields in not optimal. 



We now want to estimate the time, Tcq, required to 
build up a large scale field of equipartition field strength, 



i.e. {b') 
we have 



B^q. In units of the turnover time, r = L/un 



q/r = UrinsL/{2rikiL^eoQ) = -Rm/(eiQ), 



(17) 



where we have introduced a new coefficient ei — 2eiik\L'^ . 
Applying this to the sun we have Toq/r ^ 10''-10^ if we 
assume = 10**-10^°, Q = 10 - 100, and ei fa 2(27r)^ « 
100. At the bottom of the solar convection zone the turnover 
time is about 10 days (0.03 yr), so the time scale for building 
up a large scale field to equipartition strength is between 300 
and 3 x 10^ yr. 

We have not yet calculated models with different values 
of the magnetic Reynolds number, so we cannot properly as- 
sess the effect on the cycle period. If the cycle period scales 
in the same way as the growth time of the dynamo, then 
the helicity constraint would, even in the presence of shear, 
continue to pose a serious problem for understanding cyclic 
activity of solar-like stars. However, before making more de- 
tailed comparisons with astrophysical bodies it would be 
important to assess the importance of open boundaries, for 
example. This seems to be now one of the most important 
remaining aspects to be clarified in the theory of large scale 
dynamos; see also Blackman & Field (2000) and Kleeorin 
et al. (2000). Although initial results from simulations with 
open boundaries seem pessimistic in that respect (Branden- 
burg & Dobler 2001), the effects of open boundaries are 
likely to be more important in cases with outflows (e.g. in 
protostellar accretion discs or in active galactic nuclei). It 
should also be mentioned that large scale dynamos may op- 
erate with non-helical flows; see the recent papers by Vish- 
niac & Cho (2000) and Zheligovsky et al. 2000). This may 
relax the helicity constraint, but so far there are no simula- 
tions supporting this possibility. 



6 MEAN-FIELD INTERPRETATION 

In the absence of shear the results of the simulations could 
be modelled quite well in terms of a mean-field a^-dynamo 
with simultaneous quenching of the a-effect and the turbu- 
lent diffusivity. In this section we shall try to do the same for 
the af2-dynamo. Since the shear is strong compared with the 
inverse turnover time we can make the afi-approximation. 



i.e. we can neglect the a-effect in the equation for the gen- 
eration of the toroidal magnetic field. We also assume that 
the magnetic field varies only in the direction of the vortic- 
ity vector of the shear, which is the direction in which the 
dynamo wave travels (Yoshimura 1975). In the present case 
this is the 2-direction. Thus, the relevant equations, in terms 
of the mean vector potential A, are 



dtA^ = --SAy + -nTdlA^ 
dtAy = +aBy + rjTd^Ay 



(18) 
(19) 



where By = dzA^ and r?T = + f?t is the total (rnicroscopic 
plus turbulent) magnetic diffusivity. [In Eq. ( |l8[ ) we have 
used a particular gauge that allowed us to write the shear 
term as SAy; see Brandenburg et al. (1995) for details.] As 
in the case of the a^-dynamo, we shall assume that rjt and 
a are quenched in the same way: 



7?t = 



1 + T^BByBl^ 

t2 ^2 —2 



(20) 



where as ~ i]b is assumed, and B — B^ + By with 



Bx 



-dzAy. In the case of the a -dynamo in a periodic 



domain the two components of the magnetic field were sinu- 
soidal and phase shifted by 90° such that B was spatially 
constant. It was therefore possible to obtain the solution for 
the evolution of B in closed form. The final saturation field 
strength, Bfln, was then given by [Eq. (55) of B2001] 

(for the a^-dynamo), (21) 



Bin ^ _^ 

B^ ~' 



rjkf 



where X = aoki ~ rjToki is the kinematic growth rate of the 
-dynamo. 

In the present case of an afi-dynamo, B is no longer 
spatially constant and the solution cannot be obtained in 
closed form. We therefore resort to numerical solutions of 
Eqs (|l^)-(|2c|) using periodic boundary conditions. All the 
solutions turned out to be oscillatory with a period T, but 
the temporal structure is strongly non-harmonic; see Fig. ^ 
Note that the time dependence of Bx and By is qualitatively 
and quantitatively similar to that found in the actual sim- 
ulation (Fig. ^). The field amplitude depends on the value 
of as and agrees with that found in the simulation (Fig. ^ 
if as ft; 2. The solution also depends on the value of the 
dynamo number 



V = aokiS/{riToki 



(22) 



where t^to = rj + rjto is the kinematic value of the total tur- 
bulent magnetic diffusivity. For the model shown in Fig. ^ 
we used 2? — 10, but if D is doubled the cycle period also 
approximately doubles. Thus, D — 20 would be more repre- 
sentative for the dynamo wave at a; = (cf. Fig. ^). 

Although the present analysis is straightforward and in- 
deed quite similar to other af2-dynamos considered in the 
literature (e.g. Moffatt 1978, Krause & Radler 1980), a main 
conceptual difference is that we consider here a and rjt to 
be quenched in the same way, and that we retain the micro- 
scopic magnetic diffusion rj which is not quenched. 

We have determined the value of asBln/ B^^ and the 
cycle frequency lu = 2tt/T as a function of X/rikf for different 
values of the dynamo number "D. The results are shown in 
Fig. ^ We have checked that the different curves in Fig. |l^ 
depend only on the parameter "D, regardless of the values of 
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Figure 9. Evolution of Bx and 
dimensional mean-field model with D = 10, A = 0.015 and 
?7 = 5 X 10~*. Note that Bx has been scaled by a factor 100. 
(In this case 5 > 0, so we have plotted +Bx, and not —Bx as we 
did in Fig. | where S < 0.) 



ao, S and ?7to, provided the kinematic growth rate of the 
linearized form of Eqs (|l^) and (p^, 



A = -?7Tofci + \/ aokiS/2, 



(23) 



is kept unchanged. Note that it is now this A that enters 
in the expression X/r]ki, which we have considered as the 
control parameter for the numerical solutions displayed in 
Fig. §. 

With these preparations we can now make a detailed 
comparison with the simulation data. In the simulation the 
kinematic growth rate can be read off the first panel of Fig. S 
and turns out to be A = 0.015. Thus, with = 5 x 10~ 
and fci = 1 we have X/rjkj = 30. From Fig. we see, 
then, that uj/riTokl « 0.4. For T in the range 1000-2000 



0.006-0.003, which yields 77TO ~ 0.015-0.0075, 
we can express the 



we have to - 
respectively. 

Given the values of A and rjToki 
dynamo number as 

©/Ocrit = [A/(r?Tofc?) + l]', 

where 
find© = 

we see, then, that obB^^/B, 



(24) 

; = 2 is the critical value for dynamo action, and 
-18 for T = 1000-2000, respectively. From Fig. |l^ 



we have -Bfln ~ 0.25 and Bcq = 0.035, so B^^/B, 



60-100. In the simulations 

' '°eq ~ 50, 

and therefore q:b = 1-2, which is in agreement with the 
value inferred earlier from the field amplitude; cf. Figs ^ 
and ^ We can therefore conclude that in an aSl-dynamo, 
Q and rjt are quenched much less than in an a^-dynamo, 
where as would be around 30. If the weaker quenching for 
oscillatory afi-type dynamos is confirmed for larger values 
of the magnetic Reynolds number this would also suggest 
that the cycle period is also only weakly increased. Already 
now the cycle period is closer to the dynamical time scale 
than to the resistive one. This is best seen by comparing the 
two ratios 



ui/riToki ~ 0.4 and 



ui/riki 



6-12. 



(25) 



Note also that the values of r/To and A are very close to each 
other. This confirms again that the turbulent diffusivity is 
dynamically significant and not quenched to its microscopic 
value. 

Finally we show in Fig. hi the evolution of 
{B )/{B )max for different values of ©, and compare with 



1.0 



o 
E- 



0.1 













- - ~ ~ -a, _ 














^ 






" H 1 — 


3) = 


5 






- - 


- 2) = 


10 






A-- A- 


-A S = 


20 







10 



100 



Figure 10. Normalised saturation field strength and cycle fre- 
quency for the saturated state of a nonlinear one-dimensional 
aO-dynamo with simultaneous a- and rjt-quenchings. The diag- 
onal (dash-dotted line) in the first panel gives the result for the 
corresponding a^-dynamo (for all values of ao/^To'^l)- 
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Figure 11. Resistively dominated saturation behaviour in the 
af2-dynamo for large enough dynamo numbers (X> 20). For all 
curves we have A = 0.015 and r; = 5 X 10"''. For large values of 
"D the cycle oscillation begin to distort the curve and cause ad- 
ditional deviations from the helicity constraint (solid line), which 
is best matched for D = 20. 



the shape of the curve predicted by the helicity constraint 
of Eq. ([1^. We see that the correct shape of the helicity 
constraint is matched for D ~ 20, which corresponds to 
the value obtained for T — 2000. The fact that the helic- 
ity constrained is matched for one particular value of D is 
surprising. This suggests that in nonlinear aQ dynamo the- 
ory the dynamo number is no longer a free parameter, and 
that there is only one possible value of "D for which the he- 
licity constraint with the correct value of the microscopic 
magnetic diffusivity can be matched. 
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r CONCLUSIONS 

The present investigations have shown that the effects of 
the helicity constraint can clearly be identified, even though 
much of the field amplification results now from the shearing 
of a poloidal field. Instead of having a constraint on the 
magnetic energy in the mean field, one now has a constraint 
on the geometrical mean of the energies in the poloidal and 
toroidal mean field components. The dynamo remains time 
dependent with a typical period that is closer to a dynamical 
time scale than to a resistive one. The toroidally averaged 
field alternates in sign and shows a clear migration pattern. 

The present work has revealed that, even though the ki- 
netic helicity of the flow is near its maximum possible value, 
the poloidal field shows a great deal of 'noise', whilst the 
toroidal field does not. Power spectra of the poloidal field 
show that most of the power is in small scales, making the 
use of averages at first glance questionable. However, once 
the field is averaged over the toroidal direction the resulting 
poloidal field is governed by large scale patterns (the slope 
of the spectrum is steeper than , which is the critical 
slope for equipartition of energy between small and large 
scale fields). The presence even of a weak mean poloidal 
field is crucial for understanding the resulting large scale 
field generation in the framework of an ail-dynamo. 

Th e results of the simulations can be reproduced bv 



a mean -field gf^-dynamo where alpha-effect and turbulent 
magnetic diffusivity are quenched by the magnetic field. The 
strength of the quenching is however much weaker than 
for the corresponding a^-dynamo. The resistively limited 
growth imposed by the magnetic helicity constraint is re- 
covered for one particular value of the dynamo number. 
Whether or not the cycle period becomes catastrophically 
long in the limit of large magnetic Reynolds numbers is 
not entirely clear, because the frequency dependence shown 
in Fig. ^ seems to level off at a definite value. However, 
the value of the magnetic Reynolds number where the cy- 
cle frequency levels off shifts to larger values as the dynamo 
number is increased. If it is confirmed that large magnetic 
Reynolds numbers (based on the microscopic magnetic dif- 
fusivity) also imply large dynamo numbers (based on the 
value of the turbulent magnetic diffusivity), then the cycle 
period would probably be too long to explain the cycle pe- 
riods observed in many late type stars. On the other hand, 
using a mean-field model that obeys the magnetic helicity 
constraint we found evidence that the cycle period is con- 
trolled primarily by the dynamical time scale. 

It is important to remember that the flows considered in 
the present investigations are driven by some imposed body 
force. In astrophysical bodies the flows are driven by con- 
vection and shear. This does not directly affect the helicity 
constraint which controls the long time scales discussed here. 
However, when open boundary conditions are considered it 
may be possible that real astrophysical flows have a better 
ability to dispose small scale flelds whose magnetic helicity 
has the opposite sign of that of the large scale field. [In ex- 
ternally driven flows, open boundaries do not seem to relax 
sufficiently the constraint imposed by helicity conservation; 
see Brandenburg & Dobler (2001).] 

The driven flows considered here and in related pa- 
pers have the tremendous advantage of allowing significant 
progress to be made in understanding the simulation results 



quantitatively in terms of mean-field theory. This will be 
a much harder task for real astrophysical flows. For exam- 
ple, the helicity constraint has to our knowledge never been 
identifled in simulations of astrophysically driven flows. This 
seems to be now one of the most important tasks for future 
simulations of large scale dynamos. 
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